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We study mixtures of self-propelled and passive rod-like particles in two dimensions using Brown- 
ian dynamics simulations. The simulations demonstrate that the two species spontaneously segregate 
to generate a rich array of dynamical domain structures whose properties depend on the propulsion 
velocity, density, and composition. In addition to presenting phase diagrams as a function of the 
system parameters, we investigate the mechanisms driving segregation. We show that the difference 
in collision frequencies between self-propelled and passive rods provides a driving force for segrega- 
tion, which is amplified by the tendency of the self-propelled rods to swarm or cluster. Finally, both 
self-propelled and passive rods exhibit giant number fluctuations for sufficient propulsion velocities. 



I. INTRODUCTION 

Systems of self-propelled particles are found in nature 
on many scales, ranging from flocks of birds to colonies 
of bacteria and assemblies of cytoskeletal filaments [1- 
8]. Theoretical and experimental investigations of pro- 
pelled or otherwise active particles have established a re- 
markable array of emergent behaviors not possible in pas- 
sive systems, such as giant number fluctuations [5, 8-15], 
spontaneous flow [16, 17], and unconventional rheological 
properties [17]. Of particular relevance to our work, self- 
propelled particles have been shown to exhibit dramatic 
swarming or clustering [5, 7, 8, 10, 11, 13, 15, 18-30]. 
These unconventional behaviors arise because active par- 
ticles consume or dissipate energy to fuel internal changes 
that generate motion. Inevitably, some members of an 
active system lose this ability by breaking or dying, or 
fail to acquire it to begin with. In colonies of motile bac- 
teria, for example, there is typically a significant fraction 
of non-motile bacteria [31]. It is therefore of great inter- 
est to understand and characterize emergent behavior in 
systems composed of particles with differing motility. 

In this work, we computationally investigate the be- 
havior of an example mixed motility system - a mix- 
ture of self-propelled and passive rods interacting solely 
through excluded volume interactions. We show that, 
above threshold values of the propulsion strength and 
density of active rods, the mixed system is unstable to 
the formation of segregated domains of active and pas- 
sive rods. Depending on parameters, the domain mor- 
phology and system dynamics exhibit a rich variety of 
behavior (Fig. 1), including uniform mixing, motile clus- 
ters of active rods in a gas of passive rods, segregated 
clusters of active and passive rods, and flowing bands 
or lanes. The degree of segregation depends on the in- 
terplay between an intrinsic tendency for particles with 
different collision frequencies to demix and the emergence 
of coherent swarms or polar clusters that glide through 
a system. Our study is inspired by observations of seg- 
regation of motile and non-motile bacteria in bacterial 
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biofilms [31]. While the factors driving bacterial segre- 
gation potentially include biological mechanisms such as 
chemical signaling between bacteria, recent experiments 
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Figure 1. Emergent structures in the mixed system: (a) 
isotropic mixture, (b) swarming without coherent clusters, (c) 
coherent clusters in passive gas, (d) coherent active clusters 
with passive clustering, (e) formation of transient lanes and 
(f) a giant cluster after the lanes break down. Active rods are 
labeled red, and passive blue; black dots indicate the head of 
each active rod. L is the aspect ratio, <f> is the area fraction, 
and Pe is the Peclet Number. Movies can be found in [32]. 
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Figure 2. Representative snapshots from pure active systems. 
At low packing fraction we observe behavior consistent with 
that of previous studies [5, 7, 8, 10, 11, 13, 15, 18-30]: (a) 
small but coherent polar clusters, (b) large clusters, (c) tran- 
sient apolar swarms. New behaviors emerge at higher packing 
fractions: (d) a giant swarm exhibiting smectic and tetratic 
order, (e) lanes in a nematic prior to breakdown, (f) giant 
swarms appear after the lanes shown in (e) break down. The 
head of each rod is indicated by a black dot, and rods within 
different lanes are distinguished by color in (e). The system 
parameters are listed above each image. Movies can be found 
in [32]. 



have shown that excluded volume interactions play an 
important role in emergent behavior in bacterial suspen- 
sions, both in bulk and near boundaries [21, 33, 34], and 
theoretical work has shown that density dependent dif- 
fusion coefficients can drive pattern formation [35-37]. 
In this work, we demonstrate a generic mechanism by 
which excluded volume interactions and propulsion can 
give rise to segregation without signaling or long-range 
interactions. 



II. MODEL 

We model rods as 2D spherocylinders (rectangles with 
circular caps), represented as line segments of length t 
on a two-dimensional plane with periodic boundary con- 
ditions. We implement excluded volume interactions be- 
tween rods with a purely repulsive WCA [38] potential 



V (r) = 4/c B T 



k B T, cut off by set- 
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ting V (r) =0 for r > 2 1 / 6 , where r is the shortest dis- 
tance between the two representative line segments. This 
potential gives the rods an effective diameter a and an 
aspect ratio of L = (JLja + 1). In addition, a fraction / a 
of rods are self-propelled and hence endowed with a con- 
stant propulsive force with magnitude Fp directed along 
the rod axis. 

We do not include hydrodynamic interactions in our 
model; these can be neglected in 2D for a dense viscous 
system due to screening. We consider overdamped Brow- 
nian dynamics with the discretized equations of motion 
[39] 



r (t + St) = r (t) + H" 1 (t) ■ F s (t) St 
u(t + St) = u (t) + 7 r _1 T s x u (t) St - 
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where u is a unit vector directed along the rod axis, r is 
the position of the rod, St is the time step, F$ and T$ 
are the systematic force and torque, and Sr and Su are 
Brownian noise terms. The friction coefficients are given 
by S = 7||U (g) u + 7_l (I — u ® u), with 7^ = 27y, and 

T 2 

^7||. The Brownian noise terms are determined by 
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(St (t) Sr(t)) = 2k B TS~ 1 (t) St 

(Su (t) Su (t)} = 2k B T— (l - u (8) u) St 

7 r \ / 

The systematic force Fs = Fe + Fp includes both the ex- 
cluded volume interaction forces Fe and, for active rods, 
the uniform self propulsion force, Fp = Fpu. 

The adjustable dimensionless parameters are the as- 
pect ratio L, the propulsion force Fp = Fp (^-) , 
the area fraction 0, and the motile fraction / a . The 
Peclet number is given by Pe = LFp. We integrate 
Eqs. 1 using a second-order predictor-corrector algorithm 
[40, 41] with an adaptive time step with a maximum value 

Si = St (wr) =0-0003 [42]. 

Behavior of pure active rods: We performed simula- 
tions of TV = 6400 rods with aspect ratios L = 4, 10, 21, 
for area fractions <p G [0.05,0.8], and Peclet numbers 
Pe G [10, 120]. Initially, we set / a = 1.0 to study the be- 
havior of the purely active system. We initialized all sys- 
tems with states equilibrated without propulsion forces. 
At low Pe and (j) the system remains disordered and ho- 
mogeneous in density. Increasing Pe or <\> leads to ag- 
gregation of the rods; while long rods L = 10, 21 form 
long-lived polar clusters, the short rods L = 4 form tran- 
sient aggregates which can better be described as swarms 
(see Fig. 2). 
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Clustering of self-propelled rods has been studied ex- 
tensively [5, 7, 8, 10, 11, 13, 15, 18-30] and arises because 
collisions between self-propelled rods lead to alignment 
and correlated motion [12, 18]. Most closely related to 
our work, in a study of pure self-propelled rods Yang et 
al. [30] find that the system exhibits three qualitatively 
different cluster size distributions depending on Pe. At 
low Pe, they observe a power law distribution with an 
exponential cutoff, above a threshold Pe the distribution 
develops a plateau at large cluster sizes (as seen in the 
pioneering simulations of Peruani et al. [25]), and for 
high Pe they observe jammed immobile clusters. In com- 
parison, we will focus on low and moderate values of Pe 
and we observe only the first two forms of cluster size 
distributions (a detailed description of the cluster size 
distributions is given in the SI, Fig. SI [32]). 

We examine a larger range of packing fractions than 
has been simulated previously [25, 30] and we find that 
large packing fractions of self-propelled rods lead to some 
qualitatively new behaviors (Fig. ld-f). For example, 
long rods (L = 21) above the passive isotropic nematic 
transition <ft w 0.22 form multiple bands (Fig. le) or 
"lanes" (e.g. [43-54]). Within a band, all rods are mov- 
ing in the same direction. In contrast to human pedestri- 
ans [48] or ants [47], the rods choose their flow directions 
spontaneously. The bands appear to be transient, how- 
ever, since in most of our simulations they eventually 
break down when a small group of active rods begins 
moving perpendicular to the direction of nematic order, 
shattering the nematic order and generating configura- 
tions like that shown in Fig. If. The time before break- 
down decreases with increasing Pe, and once a system 
of lanes has broken down we never see long-lived reap- 
pearance of nematic order [55]. Another feature of our 
clusters is smectic ordering, reminiscent of behavior seen 
in vibrated granular systems [56], which we attribute to 
the high packing fraction found in these clusters. Fi- 
nally, we observe that short rods (L = 4) form swarms 
which are qualitatively different from the clusters formed 
by long rods, with significantly less polar order and more 
frequent turnover of rods. In the present work we will fo- 
cus on the behavior that emerges when passive and self- 
propelled rods are mixed; we will describe the behavior 
of pure self-propelled rods at high packing fractions in 
greater detail in a future work. 



III. SEGREGATION IN MIXED SYSTEMS 

Descriptive phenomenology: To study the mixed sys- 
tem, we performed simulations with the same ranges of 
parameters as for the purely active systems, but with 
varying motile fractions 0.1 < f a < 0.9. Except where 
mentioned otherwise, we will focus on f a = 0.5. A variety 
of emergent structures appear in this range of parameters 
(see Fig. 1). At low Pe and <p the system remains homo- 
geneous in density and composition (Fig. la). At higher 
Pe and 0, the active rods aggregate and form clusters or 
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Figure 3. Degree of segregation a over (a-c) varying aspect ra- 
tio, Peclet number, and area fraction; (d) varying Peclet num- 
ber and motile fraction. Simulations which have not reached 
steady- state are not included in the average, and parameters 
with no steady-state systems are marked in gray. Contours 
are drawn using an interpolating function. 

swarms, as in the purely active case (Fig. lb). Over a 
narrow range of Pe and (j) the diffusive timescale t& of 
the passive rods is short in comparison to the time be- 
tween collisions with clusters, and the active clusters are 
surrounded by a "gas" of passive rods, whose spatial prop- 
erties are relatively unaffected by the activity (Fig. lc). 
Upon increasing Pe or <\> collisions become fast in com- 
parison to Td and the passive rods exhibit large density 
fluctuations (Fig. Id). Essentially, the forward motion of 
active clusters forces the passive rods into dense "herds", 
leaving regions in their wake nearly vacant. Finally, for 
long rods above the isotropic-nematic transition, active 
and passive rods initially form separate lanes (Fig. le). 
In the system with / a = 0.5 each active lane is sepa- 
rated by a band of passive rods. Thus, in contrast to the 
active system where neighboring lanes flow in opposite 
directions, we did not observe a correlation between di- 
rections of motion of nearest active lanes. However, the 
active lanes do impart some extra motion to rods at the 
edges of passive lanes. As in the purely active systems, 
these lanes are unstable and ultimately break down to 
yield giant clusters of active rods (Fig. If). 

Segregation order parameter: To quantify the depen- 
dence of the degree of segregation of the active and pas- 
sive rods on the control parameters, we employ an or- 
der parameter commonly used to measure the degree to 
which two human population groups live separately from 
one another [57]. We divide the system into boxes of side 
length 6, and compute a weighted sum over boxes 

0~a (b) = Try-^ FT (~) \ f* ~ ^ ' 

2/ a (1 - /a) i Wot/ 

where rii is the total number of rods in box z, fa is the 
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motile fraction within the box, n to t = Zli n ^ an d the 
prefactor ensures that a a G [0, 1] with corresponding 
to the theoretical minimum possible segregation and 1 
corresponding to maximal segregation. It is well known 
that measurement of segregation is hindered by the mod- 
ifiable unit area problem [57] meaning that the value of 
the order parameter depends on box arrangement and 
coarse-graining length b. For example, choosing b smaller 
than the rod diameter would yield a large value of segre- 
gation even for Pe = due to the hard core repulsions, 
while choosing b comparable to the simulation box size 
would result in no apparent segregation for any value of 
Pe. We overcome this problem in two ways. First, we 
use overlapping boxes. Second, noting that the order pa- 
rameter is nonzero for Peclet number zero, we subtract 
the same quantity <Jo(b) computed for a system with the 
same parameter values except Pe = 0. The difference in 
segregation a a — o~o is nonmonotonic with respect to b (as 
shown SI Fig. S2 [32]). We report the maximum value of 
this quantity over a = max^ [<j a (b) — <Jo (b)]. Because 
a subtracts the small degree of segregation that occurs 
naturally in passive systems, its maximum possible value 
at complete segregation is approximately 0.75 (this value 
depends weakly on area fraction). 

Dependence of segregation on system parameters: As 
shown in Fig. 3, the degree of segregation a generally in- 
creases with increasing Peclet number, area fraction, and 
aspect ratio, with the large aspect ratio rods approach- 
ing perfect segregation for large Pe and (j). However, the 
degree of segregation decreases mildly for large (j) and 
small L (Fig. 3b) where it is frustrated by jamming. For 
L = 21 and low Pe, the nematic lanes did not break down 
within our simulation time (these parameter values are 
indicated by gray boxes in Fig. 3c). Finally, as shown in 
Fig. 3d, the degree of segregation is nonmonotonic with 
respect to the motile fraction / a , the maximum occurring 
roughly at the case of a symmetric mixture f a = 0.5. 
Snapshots illustrating these behaviors are found in the 
SI, Figs. S3-S6 [32]. 



IV. ANOMALOUS DENSITY FLUCTUATIONS 

While we have already noted that clustering in our 
system leads to large magnitude density fluctuations, it 
has been predicted theoretically [9, 12] and shown nu- 
merically [10, 14, 58] and experimentally [8, 11, 13, 15] 
that active systems with orientational order exhibit 'gi- 
ant' number fluctuations (GNF), meaning that the vari- 
ance grows faster than the mean as one increases the 
size of a region of interest. To understand the degree 
to which activity transfers from the self-propelled rods 
to their passive neighbors, we monitored number fluc- 
tuations of both active and passive rods. Specifically, 
we chose square regions of interest of varying size dis- 
tributed throughout the system and measured separately 
the number of active and passive particles n a (t), n p (t) as 
a function of time. We find that the scaling of num- 
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Figure 4. Density fluctuations for active and passive rods 
separately in mixed systems, and for all rods in pure active 
systems. The parameter is the exponent 5 of the scaling rela- 
tion cr n oc {n) S for the fluctuations of density in boxes which 
are a subset of the entire simulation box. The exponent is 
calculated for regions with mean number 1 < (n) < 100 to 
avoid finite size effects, shown in the bottom figure. Contours 
are drawn using an interpolating function. 

ber fluctuations with respect to the mean (n) exhibits a 
power law a oc (n) 6 for boxes with 1 < (n) < 100; den- 
sity fluctuations in larger boxes are suppressed by our 
finite system size of 6400 rods (Fig. 4). As expected, 
the active particles exhibit GNF, meaning that S > 0.5. 
Although previous experiments and simulations of self- 
propelled rods or discs found scaling exponents in the 
range 0.75-0.8 [11, 58], we find (Fig. 4) that 6 gradually 
increases with Pe and 0, and appears to saturate at about 
5 ~ 0.9, with a sharper increase for larger aspect ratios. 
Similar numbers are observed if we measure fluctuations 
for all rods, regardless of their motility state, and we ex- 
clude any systems which have not reached steady-state 
from the average. 

The fact that the density fluctuation exponent under- 
goes a gradual transition can be understood by noting 
that a theoretical estimate of S = 1 [9, 12] was derived 
for the case of an active nematic, where the active rods 
enhance density fluctuations by advecting in the direc- 
tion of alignment. We expect a different result for an ac- 
tive but isotropic system. For parameter values at which 
the measured density fluctuation exponent gradually in- 
creases with Pe and the fraction of rods incorporated 
within aligned clusters is gradually increasing. Consis- 
tent with this explanation, the increase is most gradual 
for the aspect ratio 4 rods, which would not form a ne- 
matic phase when passive. 

As shown in Fig. 4, we find that the passive rods also 
show GNF, with a similar scaling exponent 5 ~ 0.9 for 
large activities and densities. However, the value of S 
increases more slowly with activity and density than it 
does for the active rods. The value of the scaling expo- 
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nent is correlated to the amount of clustering of passive 
rods; as described above, passive rods cluster extensively 
at large values of Pe, when they collide with active rods 
frequently in comparison to the diffusive timescale. In 
contrast, for moderate self propulsion velocities collisions 
are relatively infrequent, relatively little of the activity is 
transferred from active rods to their passive neighbors, 
and thus the behavior of the passive rods resembles that 
of a pure equilibrium system (Fig. lc). 



V. ORIGINS OF SEGREGATION 

Analysis of the segregation phase diagram (Fig. 3) 
shows that multiple mechanisms control segregation in 
our system. In this section, each of these mechanisms 
and the interplay among them is discussed. 

Collision frequency disparity: First, due to their self 
propulsion, active rods experience more frequent colli- 
sions than their passive counterparts. The following sim- 
ple analysis shows that this difference in collision fre- 
quencies can provide a generic mechanism for segrega- 
tion. This is the dominant physics both at low Pe and 
low (j) and at initial times at higher packing fractions. 

Given an overdamped system of two species of parti- 
cles, the macroscopic description will be a coupled set of 
diffusion equations for the density of active and passive 
rods p a and p p respectively, 

dtPa =AiaV 2 pa + AipV 2 pp 

d t pp =£> pp V 2 pp + £> pa V 2 /V 

Linearizing the the above equations about a homoge- 
neous mixed state p a?p = p® p + £p a , P , introducing a 
Fourier transform X — J dre % ^' r X and looking for so- 
lutions of the form p = e s ( fc ^, we can identify the linear 
modes of this system as 



s± (k) 



(D aa + D pp ) fc2 



± y \A Daa ~ d pp )2 + 4D *p d p 



Clearly the + mode goes unstable whenever 



-^ap-^pa ^ -^aa-^pp- 



(2) 



This is the manifestation of the well-known phenomenon 
of cross-diffusion among species rendering the homoge- 
neous mixed state unstable. 

We can estimate the parameters in Eq. 2 as follows. (1) 
To leading order in density, the diffusion of the passive 
rods captured in D pp will correspond to the "thermal" 
diffusion coefficient of an isolated rod, D pp ~ k^T/^f = 
Do where 7 = (7N + 7j_) /2 is a mean friction coefficient. 

(2) In addition to thermal diffusion, the active rods un- 
dergo persistent motion due to their self-propulsion. This 
persistent motion is rendered diffusive at long times by 



rotational diffusion and enhances the diffusion coefficient, 
which can be estimated as [59] D aa ~ k^T/^f + Vq/Dr = 

D (1 + Pe 2 ) where we used Pe 



D D K ' 



(3) The cross diffusion terms arise from the interactions 
among active and passive rods. These can be estimated 
from the microscopic collisional dynamics by considering 
the short range repulsive interaction among the rods in 
the limit of hard particle collisions. In this limit, the 
dynamics of the density of one of the species, say the 
active rods, p a , due to collisions with the other species, 
in this case the passive rods, takes the form 



(dtp- 



a; coll 



V J 



coll 
ap * 



Here J^° n is the momentum transfer due to collisions be- 
tween particles of the two species and can be estimated 
as J^p 11 ~ ^F ap p a with F ap the mean field force acting 
on the active rods due to collisions with the passive rods. 
This force is of the form F ap (r) ~ (/<j c Ap) where / is 
the incoming flux of particles, cr c is the collision cross 
section and Ap is the momentum transfer during such a 
collision with () denoting the mean field average. The in- 
coming flux scales as / ~ P P |vi — v 2 | with the velocity 
of particle z, and the collision cross section cr c for a rod in 
two dimensions scales with its length L. Further, for mo- 
mentum conserving collisions, the momentum transfer at 

contact is of the form ^(vi — v 2 ) • k, i.e., proportional 

to the relative velocity and directed along the normal at 
the point of contact k. Finally, we must account for the 
fact that these are collisions among particles of finite size 
and therefore the relevant flux / in the evaluation of the 
mean field force should be / ~ p p (r — £) |vi — v 2 | where 
£ is the contact vector, which joins the center of mass of 
the colliding particle to the point of contact and typically 
scales as L. By assuming that the velocities of the active 
particles have a Maxwell-Boltzmann distribution of the 

form exp ^— m ^ 2/ Tjj^ ) and the velocities of the passive 

particles have a distribution exp 2^3 t ) 7 anc ^ ex P an d- 
ing the nonlocal density in gradients, the mean field force 



can be estimated to be F ap (r) ~ — L 2 (k^T - 



Vp P . 



A similar estimate can be made for the collisional contri- 
bution to the dynamics of the density of passive rods. To 
transfer this analysis from the hard particle limit to our 
overdamped con tinuous system, we equate the thermal 
velocity v t h = ^Jh^Tjm with the mean rate of displace- 
ment due to diffusion ^DqD^. Thus, the cross diffusion 



coefficients in Eq. 2 are D ap ~ (l + Pe z ) p a L' 2 and 
A- ~ ^ (1 + p e 2 ) PpL 2 respectively. 

Substituting the estimated values of these parameters 
into Eq. 2 we find that this mechanism renders the ho- 
mogeneous mixture of active and passive rods unstable 
whenever 



(1 + Pe 2 ) p aPp L 4 > 1 



(3) 
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Figure 5. Segregation parameter a versus mean cluster size. 
Top: Parameters are L = 10, / a = 0.5, varying <j> £ [0.05,0.6] 
for each Pe. The segregation a and the mean cluster size 
are strongly correlated. The points that do not fall on the 
main curve are due to the identification of clusters in the 
passive system due to random clustering. Bottom: The 
relationship between mean cluster size and segregation for 
e [0.05, 0.8] and Pe e [10, 120] changes with varying / a , but 
not as markedly with varying L. For / a close to 1.0, the re- 
lationship weakens because segregation is no longer governed 
solely by the formation of coherent clusters. 
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Figure 6. Comparison of mean cluster sizes between purely 
active and mixed systems. Holding the density of active rods 
fixed and adding passive rods enhances clustering for • f a 
less than about 0.3, but weakens it at very high area fractions 
by limiting motion. Parameters are L — 10, varying (f) £ 
[0.05,0.8] for each f a . 



It is important to note that the preceding calculation 
neglects numerical factors which are not necessarily of or- 
der 1, and can be expected to break down once clustering 
occurs since we have not accounted for polar order. How- 
ever, it does illustrate the mechanism by which dispari- 
ties in collision frequencies can lead to segregation, and 
we expect the trends predicted in Eq. 3 to be relevant. In- 
deed, the prediction that segregation occurs more readily 
with increasing Pe, packing fraction and aspect ratio is 
consistent with all of the simulation results at applicable 
parameter values. Finally, we note that analogous effects 
lead to clustering and/or lane formation in colloidal sys- 
tems for which interaction with external fields leads to 
different particle mobilities (e.g. sedimentation of parti- 
cles with different masses [44-46]) or electrophoresis of 
particles with different charges (e.g. [43]). 

Positive feedback between clustering and segregation: 
For much of the parameter space, the generic collision in- 
duced mechanism for segregation appears to be enhanced 
by the clustering that occurs in our self-propelled rod sys- 
tem. We find that the segregation behavior is closely cor- 
related with the clustering phenomenon associated with 
self-propelled rods (see Fig. 5). At the same time, seg- 
regation also induces further clustering. The alignment 
resulting from collisions of self-propelled particles can be 
thought of as an effective attraction between active par- 
ticles, limiting diffusion in the transverse direction [18]. 
Since this effect is smaller for active/passive collisions, it 
generates an additional driving force for segregation. The 
feedback between clustering and collision-induced segre- 
gation is reminiscent of the mechanism by which density- 
dependent diffusion coefficients can amplify density in- 
homogeneities to provide a route to phase separation in 
bacterial suspensions [35-37], but has a different origin. 
While the latter effect arises because bacterial motility 
is reduced in dense regions, the dominant effect driving 
segregation here is the suppression of the transverse diffu- 
sion coefficient due to alignment. The fact that diffusion 
is impeded at high densities actually inhibits segregation, 
as evidenced by the reduction in the segregation param- 
eter for high <j) in Fig. 3b. 

As shown in Fig. 5, there is a close relationship be- 
tween the extent of clustering, measured by the mean 
cluster size, and the segregation parameter a at moderate 
motility fractions / a < 0.5, with values of a from simula- 
tions at different Pe, and L collapsing onto a universal 
curve for Pe > 3. (For Pe < 1 the 'clusters' identified by 
our clustering algorithm correspond to active rods which 
transiently colocalize as in a passive system.) We show in 
the SI (Fig. S7 [32]) that this relationship can be under- 
stood by calculating the fraction of active rods which are 
in the interiors of clusters, based on the observations that 
most active rods are in well-defined active clusters and 
active clusters are 'pure', meaning that they incorporate 
only a small fraction (of order 1%) of passive rods. The 
cluster purity depends on Pe but is roughly independent 
of the overall motility fraction. 

It is important to emphasize that a tendency to cluster 
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alone is not sufficient to explain the data in Fig. 5. There 
must be a mechanism driving passive rods out of active 
clusters and vice versa in order for the geometric relation- 
ship to hold. Indeed, the relationship does break down 
for high motility fractions (see Fig. 5b), where clusters 
are less well-defined, a significant fraction of the passive 
rods are trapped within active clusters, and herds of pas- 
sive rods are transiently associated with active clusters 
by the mixing mechanism described below. 

While the preceding evidence establishes that cluster- 
ing enhances segregation, comparison between the pure 
active (/ a = 1) and mixed (/ a = 0.5) cases indicates that 
segregation also amplifies clustering at moderate packing 
fractions. This effect can be seen in Fig. 6, where mean 
cluster sizes are shown as a function of the density of 
active rods, 0/ a , for f a = l and f a = 0.5. We see that 
introducing passive rods (/ a = 0.5) yields larger clusters 
than do pure active simulations. We attribute the in- 
crease in cluster size to the fact that, once segregation 
occurs, the active rods experience an effectively higher 
density which enhances their frequency of interactions. 
Equivalently, the passive rods limit the transverse diffu- 
sivity of the active rods, which facilitates the formation of 
clusters. In contrast, at higher overall area fractions the 
introduction of passive rods decreases clustering. Passive 
rods inhibit clustering of active rods at high densities by 
hindering their propulsive motion and alignment. 

Cluster- collision induced mixing: A final mechanism 
affecting segregation is mixing due to collisions of ac- 
tive clusters. As described earlier, active clusters tend 
to drive locally dense regions of passive rods in front of 
them. During cluster collisions these passive rods can 
be trapped within clusters. Hence, frequent cluster colli- 
sions lead to a decrease in segregation and thus increased 
clustering under these conditions can inhibit segregation. 

The three mechanisms just described (collision fre- 
quency disparities, clustering, mixing) lead to a rich in- 
terplay among the motions of active and passive rods, 
segregation, and clustering that results in a nonmono- 
tonic variation of the degree of segregation with the 
motility fraction. At high motility fractions, the small 
number of passive rods in the system are incorporated 
within clusters, leading to slightly lower magnitudes of 
segregation at high / a (Fig. 3d). As the mean composi- 
tion of passive rods (1 — / a ) increases beyond this fraction, 
the generic segregation mechanism causes additional pas- 
sive rods to be expelled from the clusters. These rods 
are driven into clusters of their own as described above 



(Fig. Id). Thus decreasing the motile fractions toward 
/ a = 0.5 increases segregation. At low motile fractions, 
on the other hand, the preponderance of passive rods lim- 
its interactions between active rods and clustering, hence 
inhibiting segregation by either mechanism. 

VI. CONCLUSION 

In summary, we have characterized mixtures of self- 
propelled and passive rods. We find that activity drives 
segregation under a wide range of system densities, 
propulsion velocities and compositions. By studying a 
minimal model, we show that segregation does not re- 
quire processes such as chemotaxis, inter-organismal sig- 
naling, or hydrodynamic interactions, but rather can 
arise generically in any system in which particles ex- 
hibit variable motility. The alignment effects of the self- 
propulsion and rod geometry of the particles studied here 
enhance segregation, as revealed by the dependence of 
segregation on aspect ratio. 

Because the segregation mechanism found here arises 
from a minimal model, it can be studied in a wide va- 
riety of experimental systems. For example, our re- 
sults suggest that motile bacteria mixed with bacteria 
whose motility is genetically disabled should segregate in 
a range of densities and swimming speeds, although the 
model here would only apply for a dense two-dimensional 
system since we have not considered hydrodynamics. 
The predictions could be more directly compared to sys- 
tems of propelled synthetic rods. Shaken granular rods 
[13, 15, 56, 60-64], for example, could be examined with 
a mixture of polar and apolar rods, and rods propelled by 
catalytic reactions [65-70] could be mixed with passive 
rods on a substrate. We anticipate that understanding 
the relationship between differences in activity of sub- 
populations and their segregation within these model sys- 
tems could have implications for systems such as biofilms 
[31, 71], where these generic physical mechanisms under- 
lie complex biological interactions. 
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